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Near the sample edge, or a sharp magnetic field step the drift of two-dimensional electrons in a 
magnetic field has the form of skipping/snake orbits. We show that families of skipping/snake orbits 
of electrons injected at one point inside a 2D metal generically exhibit caustics folds, cusps and cusp 
triplets, and, in one extreme case, a section of the batterfly bifurcation. Periodic appearance of 
singularities along the ±i?-interface leads to the magneto-oscillations of nonlocal conductance in 
multi-terminal electronic devices. 
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Skiping orbits have been introduced into physics of 
metals by Niels Bohr in the early studies of diamag- 
netism pH. Skipping orbits also play a special role in 
the two-dimensional (2D) electron systems, where they 
determine the chiral current-carrying properties of the 
electron edge states [2^ important for the formation of 
the quantum Hall effect [3]. Snake orbits have been dis- 
cussed, first, in the context of the electrons propagation 
along domain walls in ferromagnetic metals [9^ J^ and, 
later, used for channeling ballistic electrons in 2D semi- 
conductors in a spatially alternating magnetic field [11]- 
[16j, B = l^^signx. In isotropic 2D metals, where all 
electrons with energies close to the Fermi level revolve 
along cyclotron circles with the same radius, R = pr/eB^ 
skipping and snake orbits have the form of consecutive 
circular segments matched by the specular reflection at 
the edge, or by the smooth continuity on the oposite 
sides of the ±5-interface, see Fig. 1. A close relation be- 
tween these two families can be established [17] by fold- 
ing a sheet of a 2D material (graphene [TF or hexagonal 
transition-metal dichalcogenide [19^^ which can sustain 
sharp bends with the radia much smaller than the elec- 
tron cyclotron radius) in a homogeneous magnetic field. 
Then, the electron path near the fold looks like a skipping 
orbit with circular segments alternating between the top 
and bottom layers, but when projected onto an unfolded 
sheet, the electron motion resembles the motion near a 
±5- interface. 

Below, we study bunching in the families of skipping 
and snake orbits of electrons injected into a 2D metal and 
singularites in the spatial distribution of electronic trajec- 
tories. Mathematically, such features originate from the 
singularities in the differentiable maps in Thom's catas- 
trophe theory [211 El] . The list of stable singularities in 
2D dynamical systems includes caustic folds (caustics) 
and cusps, whereas focus is the most famous unstable 
singularity. All of those are often encountered in ray 
optics [23] . While focusing of light is one of the widest- 
implemented physical phenomena, one often encounters 
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FIG. 1. Caustics (blue) of snake orbits (red) for: (a) R < 
Xq; (b) R > Xq. (c) Positions, v^ of cusps on the r.h.s 
(bottom) and l.h.s. (top) from the zbB- interface. Top panels 
show semicalssaically calculated interference patterns in the 
vicinity of classical singularities. 



optical caustics when observing sunlight sparkling on the 
sea, or starlight twinkling [25 . In electronics, obser- 
vations of caustics and focusing are less common. In 
2D systems, electron focusing requires designing poten- 
tial lenses |26] and electrostatic mirrors [27] by taylor- 
gating semiconductor heterostructures. Caustics and fo- 



cusing of surface-band electrons have also been observed 
in 'corals' built by the STM manipulation of atoms on 
noble metal surfaces [28 . Focusing has also been re- 
alised using a mmagnetic field. Magnetic focusing of 
electrons in the three-dimensional space has facilitated 
studies of Fermi surfaces of pure metals and semimet- 
als [29H3T] . Later, this magnetic focusing technique has 
lent its name to the nonlocal magnetotransport effect ob- 
served near the edge of a 2D electron gas in semiconduc- 
tors [IHH], where, inspirationally for this work, one family 
of caustics has been identified [5 for skipping orbits of 
electrons injected from the edge of a ballistic 2D elec- 
tron system. Here, we demonstrate that caustic bunch- 
ing is generic for snake/skipping orbits of electrons in- 
jected into a 2D metal at any distance Xq < 2R near 
the dzB- interface /sample edge. As shown in Fig. 1, for 
Xq > R > \Xq^ the networks of caustics displays a pe- 
riodic appearence of individual cusps, which split into 
cusp triplets when R > Xq. In general, such crossover 
would happen via the formation of swallowtail singulari- 
ties ^3 5 but, uniquely for a sharp field step/sample edge, 
these two regimes are separated at Xq = R hy the but- 
terfly bifurcation [23j seen as a higher-intensity unstable 
singularity of the intensity. 

For electrons isotropically injected into a 2D metal 
with an isotropic dispersion of electrons, at a point 
(— Xo,0) near the ±5-interface/edge at x = 0, their tra- 
jectories can be parametrised using the angle 6 (counted 
in the anticlockwise direction) between the initial veloci- 
ties and X-axis. In Figure 1, these trajectories are drawn 
for < ^ < 27r, with a step of 0.1. For the orbits near 
the ±5-interface/edge, these are the sequences of semi- 
circles with the coordinates r^ = (x^, ^n), 

Xn = Cn -^ Rsirnp, (n = b"^ {RsiuO - Xq) , 

yn=r]n^ R COS cp, TJn = 2n^ R^ - Q - R COS 6>, 

where n = 0, 1, 2, ... labels the number of times the tra- 
jectory arrived at the interface/edge at x = 0, angle ip 
counted from the ?/-axis allows to describe all points on a 
single segment {h^R^vnip < —(o^ foi" ^ > 0)^ and 6 = ±1 
for skipping/snake orbits. A sheet density. 
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of such trajectories can be evaluated (using a sequence 
of substitutions) , as 



p{v) = [1 — h^sigiix] 



OF 



dO 



(1) 



F=0 

F{x, y; e) = {x- C„)' + (y - Vnf - R^ ■ 



This density is singular along caustics R^ = {un-,Vn)^ 
where, simultaneously. 



dF 



0;F = 0. 



(2) 



This determines the equations for causitcs, 

/■ID ^^^ 2nCo , ^ 

^n = Cn ± R—^===, Cn = , _, - tan 6^, 



VTT4' VR^^ 



Vn=r]n±R 
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, s = sign(cos^). 



(3) 



where we choose the sign '±' and permitted range of 
using the requirement that 



The density of trajectories is most singular in the vicin- 
ity of cusps (i^^, v^), which are characterised by the con- 
dition ^^ = 0, additional to Eq. i2l. For strong mag- 
netic fields, such that Xq > R > ^Xq, the periodi- 
cally appearing cusps are illustrated in Fig. 1(a), with 
u^ = Xq and the universal curving of caustics near the 
cusps. 



.<«±K-t;„)3/2. 



(4) 



A very typical for cusps calculated semiclassically inter- 
ference pattern is shown in the top inset. The y-positions 
of the cusps, v^ are plotted in Fig. 1(c) against the ratio 
R/Xq. The top/bottom parts of Fig. 1(c) distinguish 
between the cusps appearing on the left/right from the 
magnetic field step. For skipping orbits, the latter should 
be folded onto the same half-plane. When R < ^Xq, all 
the extended caustics disappear, leaving only one limit- 
ing caustic of the closed orbits: a circle with a 2i?-radius 
centred at the source. The inset in Fig. 1(c) shows the 
limiting positions of the cusps at i? ^ ^Xq (spaced with 
the period of « 0.4Xo) before their final disappearence. 

The second regime in the formation of catastrophes 
of snake/skipping orbits is characteristic for the weak 
magnetic fields, such that R > Xq. In the latter case, 
the singularities form cusp triplets shown in Fig. 1(b). 
The semiclassically calculated interference pattern of the 
electron waves in the vicinity of the cup triplets is illus- 
trated in the top inset in Fig. 1(b), and their positions 
are shown on the r.h.s. of Fig. 1(c). Note that for 
Xq/R — > caustics in Fig. 3(b) transform into caustics 
of skipping orbits originated from a point source exactly 
at the edge of a 2D system [5 . 

Finally, the form of caustics and cusps at the transi- 
tion point between these two regimes, i.e., for Xq = R, 
is shown in Fig. 2. In this case, new higher-order singu- 
larities are formed - the batterfly bifurcation, for which 
In this case higher order singularities are formed out of 
two merging pieces of caustics. 
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Such a singularity is characterized by two additional con- 
straints, ^^ = with N = 1,2,3,4, and it represents 
a section of the A^ butterfly bifurcation [23] . According 
to the mathematical catastrophe theory [2TJ [22l |25] , such 




FIG. 2. Snake orbits (left) and skipping orbits (right) for 
electrons injected at a distance Xq = R from the zLB- 
interf ace/edge, with caustics (blue) mergin at a higher-order 
singularity: a section of the so called batterfly bifurcation 
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higher-order singularities are not stable, so that their fo- 
mation is uniquely specific to the infinitely sharp ±5- 
interface. Any weak smearing of the interface, or an 
effective gauge field created for electrons by lattice de- 
formations in, e.g., the bended region of folded graphene 
sheet [32 replaces this transition by a precursive forma- 
tion of the third-order singularity near the already exist- 
ing cusp: a 'swallowtail' catastrophe [23] consisting in the 
nucleation of a pair of cusps, followed by a gradual sepa- 
ration of the latter until the three singularities form the 
triplets shown in Fig. 1(b). All of the above results are 
also applicable to the electron skipping orbits, by folding 
caustics in Fig. 1 onto a single half-plane as shown on 
the r.h.s. of Fig. 2. 

Periodic appearence of cusps of snake/skipping orbits 
suggests that they can generate magneto-oscillation of 
non-local conductance in ballistic multi-terminal devices 
incorporating a ±5-interface (or graphene fold in a mag- 
netic field). Figure 3 shows the calculated magneto- 
oscillations of non-local current for two geometries of 
such devices. In Figure 3(a) we plot the magnetic-field- 
dependent proportion, Isd/Iq of the current /q injected 
from a point contact near the ±5-interface which reaches 
the drain on the r.h.s. of the edge. This was calculated 
by following the propagation of each of the injected elec- 
trons for the time up to lOW/vp (where W is the dis- 
tance from the source to the upper edge), which mimics 
the effect of a finite mean free path, i ~ lOW^, in the 



FIG. 3. (a) Bunching of trajectories (for Xq = R) and 
the calculated magneto-oscillations of the current Isd in the 
three-terminal geometry (for W — tt^Xq). Here, the injected 
current is registered in the drain placed on the right from the 
±B field step. The marks a^a jfi^fi' relate maxima/minima 
of Isd to the cusps reaching the upper edge of the sample 
on the right/left from the field step, as marked on Fig. 1(c). 
(b) Bunching of trajectories and magneto-oscillations of the 
source-drain current in a four-terminal device with a diB- 
interface and width W . Two top panels illustrate families 
of trajectories for the conditions corresponding to the max- 
imum/minimum [a/ (3) of IsD^ with positions of the cusps 
pointed by arrows. 



system. The oscillations in Isd are the result of singu- 
larities in the ensemble of the electron trajectories: each 
time when a cusp/focus on the right from the interface 
reaches the upper sample edge, Isd experiences a max- 
imim, and when a singularity on the left reaches the sam- 
ple edge - a minimum. Such oscillatory behavior persist 
both in the regime of individual cusps formation and the 
regime of cusp triplets. However, for the lowest mag- 
netic fields, such that R/Xq <C 1, the cusps in each of 
triplet get separated so much that one of those crosses the 
±5-interface; after that, the magneto-oscillations of Isd 
become rather irregular and loose in the amplitude. Fig- 
ure 3(b) gives an example of magneto-oscillations of the 
current Isd in a 4-terminal device incoporating a ±5- 
interface. Here, current is injected from an isotropic side 



contact at the lower edge (biased against the electrode on 
the l.h.s. at the upper edge) and registred using a drain 
contact placed at the r.h.s. at the upper edge. Similarly 
to Fig. 3(a), oscillations of Isd in Fig. 3(b) reflect pe- 
riodic appearence (on the left and right hand side from 
the ±5-interface) of cusps in the family of sequentially 
linked skipping and snake orbits. 

To summarise, we show that snake/skipping orbits of 
electrons injected at one point into in a 2D metal (at 
the distance Xq from the ±5-interface/edge) generically 
display caustic bunching and formation of intense local 
singularities - cusps, with two characteristic regimes of 
cusp formation: (i) periodic appearence of individual 
cusps (for 2R > Xq > R) and (ii) cusp triplets (for 
R > Xq). Singularities in the distribution of trajecto- 
ries, which are most intense when R = Xq, can lead to 
the magneto-oscillations in the non-local conductance of 
multi-terminal devices made, e.g., of a bi- folded graphene 
flake. Alternatively, one can employ near-field optics 
source to generate electron-hole pairs in the heterostruc- 
ture, with electrons placed at the energy close to the 
Fermi level, and, then, detecting the presence of singu- 
larities by measuring magnetic-field and source-position 
dependences of a voltage drop between a fixed point con- 
tact and a massive contact placed further up along the 
edge. 
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